System  Reliability  Simulation  and 
Optimization  by  Component  Reliability 

Allocation 


Zissimos  P.  Mourelatos 

Professor  and  Head 
Mechanical  Engineering  Department 
Oakland  University 
Rochester  Ml  48309 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

23  NOV  2011  Briefing  Charts 

3.  DATES  COVERED 

23-11-2011  to  23-11-2011 

4.  TITLE  AND  SUBTITLE 

SYSTEM  RELIABILITY  SIMULATION  AND  OPTIMIZATION  BY 
COMPONENT  RELIABILITY  ALLOCATION 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Zissimos  Mourelatos 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

Oakland  University, Mechanical  Engineering 

Department, Rochester, MI, 48309 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

;  #22117 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  TARDEC,  6501  E.ll  Mile  Rd,  Warren,  MI,  48397-5000 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

TARDEC 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

#22117 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

SimBRS  Annual  Review  Meeting  2011 

14.  ABSTRACT 

NA 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ABSTRACT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  3S 

unclassified  unclassified  unclassified  Report  (SAR) 

12 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Reliability  Allocation  and  Cost  Tradeoffs 

using  Simulation 

Reliability  allocation  for  a  system  is  the  process 
of  assigning  component  reliabilities  to  achieve  a 
system  level  reliability. 

Reliability  comes  at  a  cost. 

The  objective  of  this  study  is  to  simulate  an 
engineering  system  using: 

-  Newer  simulation  methods 

-  Developing  better  failure  detection  techniques 

-  Perform  multiobjective  optimization 

-  Allow  for  post-optimal  tradeoffs 
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Four  Joint  Robot  Example 

The  system  we  use  is  a  four-joint  robot  manipulator 
(Talon  IV  Engineer  arm)  consisting  of  four  motors  and 
eight  sensors  (two  for  each  joint).  One  motor  has 
redundancy  so  we  have  five  motors  effectively.  Here  is 
kinematic  diagram. 
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Four-Joint  Robot  Example  -  Fault  Tree 


robot  fails 
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Adjacency  matrix  for  the  system 
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Determining  failure  of  the  system 


Packet  travels  upstream  one  step  at  a  time.  v0  is  the  vector 
that  counts  the  copies  of  the  packet  at  each  component  in  the 
beginning 

Vi  =  Rv  o 


Where  v.,  counts  the  copies  of  the  packet  at  each  component 
after  step  1 .  To  see  if  will  eventually  reach  Cso,  multiply  v0 
successively  with  R,  at  most  n+ 1  times.  Therefore  after  p 
steps  vp  will  be  given  by. 


v  =RP 

p 


0 


Since  we  need  to  find  if  Cso  ever  gets  the  packet,  we  need  to 
sum  all  v/s  and  see  if  the  first  element  of  the  sum  is  non-zero. 
Let  this  sum  be  s: 


n+ 1 


i= 1 
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Determining  failure  of  the 
system  (contd.) 

In  case  some  of  the  components  fail,  we  can  incorporate 
this  into  the  adjacency  matrix. 

R-Rx  diag( c) 

The  elements  of  the  vector  c  are  all  one  except  those 
corresponding  to  components  that  have  failed.  These 
entries  are  zero.  The  sum  becomes: 

n+ 1 

i= 1 

RHS  is  a  geometric  progression.  We  have, 

s  =  (i  -  R')-1  v 
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Implementing  the  method 

We  introduce  a  packet  at  CSi,  this  can  be  represented  by  a 
column  vector,  all  of  whose  elements  are  0  except  the  last  one 
which  is  1: 

v0=[0  0  0  0  •••  if 

The  sum  s,  is  found  using  the  expression  shown  before. 

The  first  element  of  s,  corresponding  to  Cso  is  an  indicator  of 
failure: 

•  IF  the  first  element  is  0,  the  system  has  failed 

•  IF  the  first  element  is  1,  the  system  survives  but  has  been 
rendered  serial. 

The  number  can  also  be  interpreted  as  the  redundancy 
remaining  in  the  system.  For  example,  If  the  first  element  is  5, 
there  are  5  paths  remaining  from  CSi  to  Cso. 
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Parameters  used  for  simulation 


Baseline 
MTTF  (Hrs) 

Baseline 
Cost  ($) 

cov 

B-factor 

Motor  1 

31519 

150 

0.5 

4 

Sensor 1A 

4845 

50 

0.5 

4 

Sensor  IB 

4845 

50 

0.5 

4 

Motor  2A 

3476.62 

61.62 

0.5 

4 

Motor  2B 

3476.62 

61.62 

0.5 

4 

Sensor  2A 

4845 

50 

0.5 

4 

Sensor  2B 

4845 

50 

0.5 

4 

Motor  3 

31519 

150 

0.5 

4 

Sensor  3A 

4845 

50 

0.5 

4 

Sensor  3B 

4845 

50 

0.5 

4 

Motor  4 

31519 

150 

0.5 

4 

Sensor  4A 

4845 

50 

0.5 

4 

Sensor  4B 

4845 

50 

0.5 

4 

*Component  and  system  times-to-failure  are  assumed  to  follow  a 
Beta  distribution  because  of  its  flexibility. 
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System  Reliability  @  5000  Hrs 


Reliability-Cost  Pareto  Front 
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Designs  Associated  with  Points  1  and  2 
(Only  shaded  components  are  designed) 


MTTF  (Hrs) 

Cost  ($) 

Motor  1 

25951.9 

125.7 

Sensor 1A 

4845 

50 

Sensor  IB 

4845 

50 

Motor  2A 

4143.1 

74.6 

Motor  2B 

6203.5 

135 

Sensor  2A 

4845 

50 

Sensor  2B 

4845 

50 

Motor  3 

16128.9 

92.1 

Sensor  3A 

4845 

50 

Sensor  3B 

4845 

50 

Motor  4 

15759.5 

91 

Sensor  4A 

4845 

50 

Sensor  4B 

4845 

50 

COST 

$918.40 

RELIABILITY 

0.9125 

MTTF  (Hrs) 

Cost  ($) 

Motor  1 

18625 

99.6 

Sensor 1A 

4845 

50 

Sensor  IB 

4845 

50 

Motor  2A 

2011.3 

40.4 

Motor  2B 

2873.4 

51.8 

Sensor  2A 

4845 

50 

Sensor  2B 

4845 

50 

Motor  3 

16225.5 

92.3 

Sensor  3A 

4845 

50 

Sensor  3B 

4845 

50 

Motor  4 

17531.3 

96.2 

Sensor  4A 

4845 

50 

Sensor  4B 

4845 

50 

COST 

$780.30 

RELIABILITY 

0.6713 
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Conclusions 

Demonstrated  simulation  of  reliability  of  a  four 
joint  robot 

Developed  method  for  failure  identification  using 
partial  system  failures 

Performed  multiobjective  optimization  using 
state  of  the  art,  NSGA-II  which  gave  a  Pareto 
front 

Future  work 

-  How  to  determine  the  knee  region  (best  reliability  with 
minimal  increase  in  cost)  ? 

•  Should  use  a  utility  function  -  it  will  automatically  take  into 
account  tradeoffs 
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